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Spectral weight contributions of many-particle bound states and continuum 
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Cluster expansion methods are developed for calculating the spectral weight contributions of mul- 
tiparticle excitations - continuum and bound states - to high orders. A complete 11th order calcu- 
lation is carried out for the alternating Heisenberg chain. For A = 0.27, relevant to the material 
Cii(iV 03)2.2.51)20, we present detailed spectral weights for the two-triplet continuum and all bound 
states. We also examine variation of the relative weights of one and two-particle states with bond 
alternation from the dimerized to the uniform chain limit. 



PACS numbers: 75.40. Gb, 75.10.Jm, 75.50.Ee 

In recent years there has been a growing interest in 
understanding quantitatively the single and multiparti- 
cle excitation spectra in quantum spin systems. On the 
one hand, this understanding is important in assessing 
the extent of a continuum due to conventional many- 
particle excitations as compared to more exotic scenar- 
ios On the other hand, the increased frequency and 
wave- vector resolution of inelastic neutron scattering ex- 
periments means that such spectra are being observed in 
real materials and need quantitative theoretical support 
§• 

Controlled and systematic calculation of the spin dy- 
namics of quantum spin models remains a challenging 
computational task. Despite much progress in develop- 
ing computational methods the multiparticle excitations 
remain poorly understood. 

Here, we develop a general linked-cluster formalism to 
calculate the single-particle and multi-particle contribu- 
tions to the dynamical structure factor by means of high- 
order series expansions. We apply the method to the 
alternating Heisenberg chain (AHC) model, where ex- 
pansions are made around the strong coupling limit of 
decoupled spin dimers. 

We follow the formalism of Tennant et al. The 
inelastic neutron scattering cross-section || at temper- 
ature T — is proportional to the neutron scattering 
"structure factor" S^k.w) 

S a0 (k,u;) = J^J2j * exp[i(wt + k ■ (r, - r,)] 



tyo|S?(t)Sf(0)|^o> 



(1) 



where k is the wavevector transfer, u> is the energy, N 
is the number of scattering centres, a, (3 — x,y,z are 
Cartesian coordinates, i, j label sites of the system, and 
|^o) is the ground state of the Hamiltonian. This is just 
the space and time Fourier transform of the spin-spin 
correlation function. For AHC, the structure factor has 
only one independent spin component, and henceforth we 
concentrate our attention on S '"(k, u>). 

Integrating equation (Q) over energy, we get the inte- 
grated structure factor, i. e. the Fourier transform of the 



spin-spin correlation function at t = 0: 

S~ + (k) = J] C-+(r) exp[ik • r] 



(2) 



Inserting a complete set of momentum eigenstates \ip n ) 
of H between the spin operators in equation ([!]) , and us- 
ing translation invariance, we can express the spin struc- 
ture factor as a sum over "exclusive" structure factors 



s-+(k, W ) = ys- 



N, 



"(k,w), 



(3) 



S~ + (k,u;) = 1 ±5(u J -E n + E ) 



^^„(k)|5i(0)|^o)exp[ikT i 



(4) 



where the sum runs over all sites i* in the unit cell, and 
N c is the number of unit cells on the lattice. Barnes et 
al. [0 define the "reduced exclusive structure factor" as 



S-+(k) = N c \ £>„(k)|$+|Vto) ex p[* k • n 



(5) 



Each exclusive structure factor S""(q, ui) gives the inten- 
sity of scattering from |^>o) to a specific triplet excited 
state |V>«(k)> §. 

We turn now to a discussion of algorithms for the cal- 
culation of exclusive structure factors within perturba- 
tion theory. Efficient linked cluster expansion methods 
have long been known p] for calculating bulk proper- 
ties of a quantum lattice system. Similar methods for 
the calculation of 1-particle spectra were developed by 
Gelfand ||, and were recently extended to 2-particle 
spectra by Trebst et al Series calculations of ex- 

clusive 1-particle structure factors or spectral weights 
appear in several places (e.g. P JlO| , pT| ); and some low- 
order calculations for 2-particle states have recently been 
made P,p^|j2[ . Knetter et al j^] have used an alternative 
approach based on 'continuous unitary transformations' 
which is also capable of giving bound state energy spectra 
and structure factors to high order and in great detail. 
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Let us suppose that the Hamiltonian can be decom- 
posed 



H = Hn 



XV 



(6) 



where Hq is the unperturbed Hamiltonian and V is to 
be treated as a perturbation. We aim to expand the 
multiparticle dispersion relations and structure factors 
in powers of the parameter A. For illustrative purposes, 
we shall use the language of the AHC, but the formalism 
can be applied more generally. 

At zeroth order, the 'single-particle' excitations consist 
of triplets on a single dimer which can be labelled \m, a), 
where m labels the position of the excited dimer, and 
a labels the angular momentum eigenstate (i.e. S 2 ,S Z ). 
The reduced exclusive structure factor is 

S ~ + ( k ) = | E (m,a\St*\0)exp[ik( Xi * ~ x m )]^ (7) 

m.z* ,a 

It is easy to see that the matrix element (to, a|>Sj |0) = 
Sf a (m,i) must obey a simple 'cluster addition' property, 
such that the elements Sf a (m,i) admit a linked cluster- 
expansion 



Sf a {m,<, 



79(m,i) 



'la 



(to, i) 



(8) 



where the sum over 7 denotes a sum over all connected 
clusters which are "rooted to" , or contain, the positions 
of dimer m and site i. Correspondingly, the perturba- 
tion series expansion for S^ a (m, i) could be formulated in 
terms of a diagrammatic expansion where only connected 
diagrams contribute, although we will not elaborate on 
this approach here. 

An efficient linked cluster algorithm for calculation of 
the structure factors can now be formulated, following 
Trebst et al. 

i) Generate a list of connected clusters 7 appropriate 
to the problem at hand (in the present case, they will 
simply consist of chains of dimers of different lengths); 

ii) For each cluster 7, construct matrices for the Hamil- 
tonian H and spin operators Sf in the basis of singlet 
and triplet dimer states corresponding to H ; 

hi) 'Block diagonalize' the Hamiltonian by an orthogo- 
nal transformation as outlined by Trebst et al @^], con- 
structed order-by-order in perturbation theory so that 
the 1-particle states sit in a block by themselves; simul- 
taneously, transform the matrices for the spin operators 



err, 



(9) 



iv) Subtract all sub-cluster contributions; 

v) Insert in equation (^), and then equation (^), to 
build up the exclusive structure factors order-by-order in 
perturbation theory. 



Since two triplets can combine to give total spin 1, 
there will also be non-zero spin structure factors for 2- 
particle states. The 2-particle states can be labelled ac- 
cording to their unperturbed counterparts in position 
space |to, n; a, 0) where to, n label the two dimer po- 
sitions, and a,P the corresponding angular momentum 
states. Then the exclusive matrix element is 

(^■|5+|0> =^/^(TO,n)(TO,n;a,/3|S+|0) (10) 



where the 2-particle wavefunctions /£g(TO, n) can be com- 
puted by the block diagonalization procedure of Trebst 
et al. 0,^]. The 2-particle matrix elements in position 
space can then be found by a linked cluster algorithm 
essentially the same as that for the 1-particle states. 

We have used this method to investigate the spectral 
weights for the alternating Heisenberg chain Hamilto- 
nian, described by the following Hamiltonian 



H 



(11) 



where the Si are spin-i operators at site i, and A is the 
alternating dimerization. 

There is a considerable literature on this model, which 
has been reviewed recently by Barnes et al. j|. At A = 0, 
the system consists of a chain of decoupled dimers, and 
in the ground state each dimer is in a singlet state, while 
excited states are made up from the three triplet excited 
states on each dimer, with a finite energy gap. At A = 1, 
we regain the uniform Heisenberg chain, which is gapless. 
Numerical studies of the model include series expansions 
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11 1, and exact diagonalizations for finite lattices 



The 2-magnon bound states were previously studied 
by Uhrig and Schulz JTt]] using an RPA approach. They 
found a singlet bound state below the 2-particle con- 
tinuum for all momenta k and over the whole range of 
A < 1. They also predicted a triplet bound state and a 
quintuplet antibound state near k = 7r/2 for A not too 
large. These conclusions were supported in later studies 
@|o|Jl§. Trebst et al. @@ found in a high-order series 
expansion study that in fact there are two S = (Si and 
S 2 ) and two S = 1 bound states (Ti and T 2 ), together 
with two 5 = 2 antibound states (Qi and Q 2 ), if one goes 
sufficiently close to k = tt/2. Hence the model displays 
some interesting multi-particle dynamics which can be 
explored both theoretically and in experiments. Neutron 
scattering experiments, however, will only be sensitive to 
the triplet bound states. 

Barnes et al. and Tennant et al. |2j have shown how 
to calculate exclusive structure factors for the 2-particle 
states by low-order series expansions, and have made a 
comparison with experimental data |^l],^| for the copper 
nitrate material, Cu(A0 3 ) 2 .2.5D20- 
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FIG. 1. The integrated structure factor S, and some rela- 
tive spectral weights at A = 0.27, as functions of wavevector 
k, where d is the distance between dimers. 



Here we extend the series expansions to high orders. 
Series have been computed up to order A 11 for the inte- 
grated structure factor S(k), and the exclusive structure 
factors for the 1-particle triplet state Si p (k), the two 2- 
particle triplet bound states, St x and St 2 , the 2-particle 
continuum S2 pc , and finally the total 2-particle structure 
factor (the sum rule) S 2p - The 1-particle structure factor 
has been computed to order A 3 by Barnes et al. (J, but 
our series disagree with their results from second order. 

The integrated structure factor for A = 0.27, relevant 
to the material Cu(NOs) 2 .2.5D 2 0, is shown as a func- 
tion of wavector k in Fig. ED, along with some relative 
spectral weights. It can be seen that the structure fac- 
tor peaks strongly at kd = 2-7T, and is dominated by the 
1-particle state, which carries at least 96% of the total 
weight. The relative weight for more than two particles 
(i.e. (5 - Sip - S 2p )/S ) is less than 0.0069%. 

The spectral weights for 2-particle triplet bound states 
Ti and T 2 are given in Fig. ||. We know from previous 
studies that bound states exist only in a limited range of 
momenta near kd = n, and the spectral weight reflects 
this fact. The maximum relative spectral weight for Ti 
is about 1.2% at kd = ir, while the spectral weight for T 2 
has two peaks, with a zreo in between them, and has a 
very small maximum relative weight about 0.023%. The 
two-peaked structure can be traced back to the bound- 
state wave function for A — > 0, where the two triplets are 
separated by an odd number of dimers ^cj . 

Figure || shows the structure factor for the 2-particle 
continuum, also at A = 0.27, sliced at various momentum 
intervals. Also shown as a solid curve is the dispersion 
relation for the triplet bound state T\ . It can be seen that 
where the bound state merges with the continuum, a very 
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FIG. 2. The relative spectral weight of the 2-particle triplet 
bound states Ti and T 2 at A = 0.27. 
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FIG. 3. The structure factor (shifted by n/100) for the 
2-particle continuum versus energy u> at A = 0.27 and 
kd = 27m/20, n = 0, 1, 2, • • ■ , 20. 



sharp peak develops at the lower edge of the continuum. 

Figure |] shows first, the position of the peak in the 
continuum structure factor relative to the lower edge, 
Lu p — u>l, and secondly the peak value of the structure 
factor S2 pc {oJp,k), as functions of k, calculated with fi- 
nite lattices of 600, 1200 and 2400 sites in the Fourier 
transform [g) . It can be seen that at the threshold points 
where T± and T 2 merge with the continuum, (u> p — lul) 
goes quadratically to zero, and the peak value S^p^Wp, k) 
shows spikes which appear to be diverging to infinity for 
the bulk system. 

We have also studied the behaviour of the spectral 
weight as a function of A. Figure |E| shows, for kd = it, 
the integrated structure factor S, the relative spectral 
weight of the 1-particle triplet state, the total 2-particle 
spectral weight, and those for the individual 2-particle 
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kd/7T 

FIG. 4. The peak value of the 2-particle continuum struc- 
ture factor S% (u) p , k)' an d t fle difference (uj v — lul) between 
the position of the peak and the lower edge of the continuum, 
graphed as functions of k for finite lattices of 600 (dotted line) , 
1200 (dashed line) and 2400 (solid line) sites at A = 0.27. 



triplet bound states Ti and T 2 and continuum. It can 
be seen that the relative spectral weight of the 1-particle 
triplet state approaches zero as A — > 1, as expected, and 
the spectral weight for T 2 vanishes before A = 1, while 
that for T\ vanishes at A = 1. The 2-particle continuum 
has about 54% weight at A = 1 , implying a large overlap 
between the two-triplet and the two-spinon descriptions 
p2| in the uniform limit. 

These results show that detailed and accurate results 
for structure factors and spectral weights can be obtained 
using high-order series expansions, as was also demon- 
strated for the spin-ladder model by Knetter et al. [jl2j. 
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